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PREFACE 



This publication describes a number of 
facilities which provide instructions to 
improve performance in certain mathemat- 
ical computations. The description is 
in three parts. The first part 
describes the multiply-and-add facility. 
The second part covers the square-root 
facility. The third part describes 
facilities which evaluate certain math- 
ematical functions; they are the 
arctangent, common- logarithm, exponen- 
tial, natural -logarithm, raise-to-power, 
and sine-cosine facilities. 

The instructions are valid in any archi- 
tectural mode (System/370, ECPStVSE, or 
370-XA) available on a model equipped 
with the corresponding mathematical- 
assist facility. The reader should be 
familiar with the IBM System/ 370 Princi- 
ples of Operation, GA22-7000, the IBM 
4300 Processors Principles of Operation 
for ECPS-.VSE Mode, GA22-7070, or the IBM 
System/370 Extended Architecture Princi- 
ples of Operation, SA22-7085, as appro- 
priate, and particularly with Chapter 9, 
"Floating-Point Instructions," of any of 
those publications. 



The facilities discussed in this publi- 
cation are not available on every model. 
At the time of publication, they are 
available only on some models of the IBM 
4341, 4361, and 4381 Processors, but all 
the facilities are not provided on every 
one of those models. Publication does 
not imply any intention by IBM to 
provide the facilities on models other 
than those for which they are announced. 
For current information concerning the 
availability of the facilities on any 
specific model, refer to the latest 
edition of the functional character- 
istics publication for the model. 



Terminology 

As used in this publication, a floating- 
point scalar is a single floating-point 
number. A floating-point vector is a 
linearly ordered collection of floating- 
point numbers, each number being an 
element of the vector. Consequently, a 
vector consists of a set of elements, 
each of which is a scalar. 
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MULTIPLY-AND-ADD FACILITY 



OVERVIEW 

The MULTIPLY AND ADD instruction per- 
forms a combination of vector multipli- 
cation and addition operations which may 
replace the inner loop of common matrix 
computations. Its function may be 
described as: 

A = (B * S) + C 

where * indicates multiplication. B is 
a vector that is multiplied by the 
scalar S. The product is added to the 
vector C, and the sum replaces vector A. 

The three vectors (A, B, and C) are in 
storage. Each consists of one or more 
floating-point numbers called the ele- 
ments of the vector. Each vector con- 
tains the same number of elements . 
Scalar S is a floating-point number pre- 
viously loaded into floating-point reg- 
ister 0. The floating-point numbers are 
all in normalized form and in the long 
format of 64 bits, except that vector C 
may contain unnormalized elements. 



separation for vector C is always the 
same as for vector A, but it may differ 
from the element separation for vector 
B. 

For example, consider an N-by-N matrix 
that is stored in column order, the con- 
vention used for IBM System/ 3 70 FORTRAN 
programs. The elements of a column 
vector are contiguous, and the column 
vector has an element separation of 8. 
The elements of a row vector, however, 
are not contiguous, and a row vector has 
an element separation of 8N. The vector 
of elements along the major diagonal of 
the matrix has an element separation of 
8(N + 1). All three vector types 
contain N vector elements. 

Vectors A, B, and C all may be dif- 
ferent, if their storage locations do 
not overlap; or any two or all three may 
coincide. If either of the two source- 
operand vectors partially overlaps the 
result vector in storage, the result is 
undefined. 



Vectors in Storage 

All elements of the vectors must be 
located in storage on doubleword bounda- 
ries, so that their addresses are multi- 
ples of 8. Successive elements of a 
vector are uniformly spaced; they may be 
contiguous (in successive doublewords) 
or separated by other data. The same 
number of elements are processed in each 
vector. 

The increment in bytes from the address 
of one vector element to the next is 
called the element separation. For con- 
tiguous elements, the element separation 
is 8. If the elements are not contig- 
uous, the element separation must be a 
constant multiple of 8. The element 



Instruction Execution 

The MULTIPLY AND ADD instruction per- 
forms a sequence of operations that is 
essentially equivalent to the execution 
of the following floating-point 
instructions on each set of corre- 
sponding vector elements: 



LOAD (LD) 
MULTIPLY (MDR) 

ADD NORMALIZED (AD) 
STORE (STD) 



Load element of B 
Multiply by 
scalar S 

Add element of C 
Store result in A 



Arithmetically, the result is the same 
as if those instructions were embedded 
in a simple loop that also included 
instructions to increment the storage 
addresses from one vector element to the 



Multiply-and-Add Facility 1 



next, and an instruction to branch back 
until all elements have been processed. 

The MULTIPLY AND ADD instruction differs 
from such a loop in that the recognition 
of an exceptional arithmetic condition 
(exponent overflow, exponent underflow, 
or significance loss) does not cause a 
program exception to be recognized and 
an interruption to occur, even though 
the program mask in the PSW may permit 
the interruption. Instead, the occur- 
rence of such a condition causes 
instruction execution to be completed 
and a nonzero condition code to be set; 
the result for the current set of ele- 
ments is not stored, and no more ele- 
ments are processed. 

The MULTIPLY AND ADD instruction also 
differs when the instruction encounters 
an unnormalized multiplication operand 
or a vector element that is not located 
on a doubleword boundary in storage; 
execution is completed without storing 
the result element, and a nonzero condi- 
tion code is set. 

The MULTIPLY AND ADD instruction should 
be followed by a BRANCH ON CONDITION 
instruction to test for a zero condition 
code, indicating normal execution. 
Detection of a nonzero condition code 
may be used to cause execution of a 
series of floating-point instructions, 
other than MULTIPLY AND ADD, which refer 
to the same storage addresses and 
reprocess the current vector elements. 
If an exception is then encountered, an 
interruption occurs, which allows the 
usual action to be taken. No new 
exception-handling programs are 
required. See Programming Note 2 at the 
end of the instruction description for 
an example. 

The MULTIPLY AND ADD instruction may be 
interrupted during execution for other 
causes, such as access exceptions, and 
I/O or external interruptions. As with 
the COMPARE LOGICAL LONG (CLCL) and MOVE 
LONG (MVCL) instructions, the instruc- 
tion, when reexecuted, resumes at the 
point of interruption. 



INSTRUCTION DESCRIPTION 

The inclusion of an instruction mnemonic 
in the description does not necessarily 
imply that the mnemonic is recognized by 
current assembler programs. A pro- 
grammer wishing to use this instruction 
in assembler- language programs with an 
assembler which does not recognize the 
mnemonic may use an appropriate macro to 
assist with the translation of the 
instruction into machine language. 



MULTIPLY AND ADD Instruction 

MADS Ri,R 2 [RRE] 



| *B243' I////////I Ri I R 2 I 

l I i i I 



16 



24 



28 31 



The MULTIPLY AND ADD instruction per- 
forms the vector multiplication and 
addition operations: 

A = (B * S) + C 

where A, B, and C are three vector oper- 
ands, S is a scalar, and * indicates 
multiplication. The vector operands are 
in storage, where they must be aligned 
on doubleword boundaries. The scalar 
operand and all vector elements are 
floating-point numbers in the long 
format. The scalar and the elements of 
operand B must be in normalized form; 
the elements of operand C may be in nor- 
malized or unnormalized form. The ele- 
ments generated for operand A are 
normalized. 

The scalar operand is in floating-point 
register 0. The vector operands in 

storage are specified by the contents Of 

as many as six general registers. Three 
of the general registers have a fixed 
assignment, and any others are desig- 
nated by the R a and R 2 fields of the 
instruction. 
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General register 1 contains a 32-bit 
unsigned binary integer, which repres- 
ents the number of elements in each 
vector. General registers 2 and 3 
specify the address of the first element 
of operands A and B, respectively. 

The Rj field, if nonzero, designates a 
pair of general registers, called the 
even Rj register (numbered Rj ) and the 
odd R : register (numbered Ri+1) . The 
even Rx register contains the element 
separation for operands A and C, and the 
odd Rj register contains the element 
separation for operand B. The element 
separation for a vector is the increment 
in bytes from the address of one vector 
element in storage to the address of the 
next element. If the Ri field is zero, 
however, the field does not designate 
general registers and 1. Instead, the 
elements of all three vector operands 
are specified as contiguous in storage, 
and element separations of 8 are 
implied. The R t field must be zero or 
contain an even number; otherwise, a 
specification exception is recognized. 

The R 2 field, if nonzero, designates a 
general register, called the R 2 reg- 
ister, which contains the address of 
operand C. If the R 2 field is zero, 



however, the field does not designate 
general register 0; instead, the address 
of operand C is the same as for operand 
A, as specified by general register 2. 

The number of general-register bit posi- 
tions that are used to specify an 
address or element separation depends on 
the mode of operation. The System/370 
and ECPS:VSE architectural modes have 
only 24-bit addressing, whereas the 
370-XA architectural mode offers the 
choice of 24-bit or 31-bit addressing. 
When 24-bit addressing is in effect, bit 
positions 8-31 are used as a 24-bit 
address or element separation, and the 
contents of bits 0-7 are ignored. When 
31-bit addressing is in effect, bit 
positions 1-31 are used as a 31-bit 
address or element separation, and the 
contents of bit are ignored. 

Figures 1 and 2 illustrate the contents 
of the general registers for the two 
types of addressing, where GR1, GR2, and 
GR3 represent general registers 1, 2, 
and 3, respectively. Figure 1 shows the 
register layout for 24-bit addressing 
(System/370, ECPSrVSE, or 370-XA). 
Figure 2 shows the layout for 31 -bit 
addressing (370-XA only). 
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Figure 1. General-Register Assignment 
with 24 -Bit Addressing 
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Figure 2. General-Register Assignment 

with 31 -Bit Addressing 



4 IBM System/370 Mathematical Assists 



Figure 3 summarizes the sources of the 
operand addresses and the element sepa- 
rations, according to whether the R : and 
R2 fields are zero or nonzero. 
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R, R; 
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=0|?*0 
#0 |=0 
#01 #0 



Source of 
Address 



A I B 



(2)|(3)|(2) 
(2)|(3)|(R 2 ) 
(2)|(3)|(2) 
(2)|(3)|(R 2 ) 



Element 
Separation 



"T 



A and C| B 



8 

8 

(Ri) 
(Ri) 
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Explanation : 

(R) Contents of general 
register R 



Figure 3. Operand Addresses and Element 
Separations 



Execution of the instruction begins with 
three tests performed in the following 
order: First, if general register 1 
contains zero, condition code is set. 
Next, if the fraction of the operand in 
floating-point register is nonzero but 
unnormalized (the leftmost hexadecimal 
digit is zero), condition code 2 is set. 
Finally, bits 29-31 of all general reg- 
isters that contain vector-operand 
addresses and element separations are 
tested for zeros; if any of those bits 
is one, condition code 3 is set. 

If any one of the preceding tests sets 
the condition code, the remaining tests 
are omitted, instruction execution is 
completed, and register and storage con- 
tents remain unchanged. Otherwise, the 
operation proceeds by repeating the fol- 
lowing steps until instruction execution 
is completed or interrupted: 

1. If the fraction of the element at 
the address of operand B is nonzero 
but unnormalized (the leftmost 
hexadecimal digit is zero) , condi- 
tion code 2 is set, and instruction 



execution is completed. Otherwise, 
instruction execution continues . 

If instruction execution continues, 
the element of operand B is multi- 
plied by the contents of floating- 
point register 0. If either number 
to be multiplied has a zero frac- 
tion, the product is set to a true 
zero. The product and the element 
at the address of operand C are then 
added. Register and storage con- 
tents remain unchanged during this 
step. 

If an exponent -overflow condition 
occurs during either the multipli- 
cation or the addition, it is not 
treated as a program exception. 
Instead, condition code 1 is set, 
and instruction execution is com- 
pleted. 

If an exponent -underflow condition 
occurs during either the multipli- 
cation or the addition, or if the 
addition produces a zero result 
fraction, the condition is not 
treated as a program exception. 
Instead, instruction execution 
depends on the settings of the 
exponent -underflow and significance 
masks in the PSW: 

• When the mask that corresponds 
to the recognized condition is 
one, condition code 1 is set, 
and instruction execution is 
completed. 

• When the exponent -under flow mask 
is zero, exponent underflow 
during the multiplication causes 
a true zero, instead of the 
product, to be added to the 
operand-C element. Exponent 
underflow during the addition 
causes a true zero to replace 
the result of the addition, and 
instruction execution continues. 

• When the significance mask is 
zero and the addition produces a 
zero result fraction, a true 
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zero replaces the result, and 
instruction execution continues. 

Exponent overflow or exponent 
underflow during the multiplication 
is recognized even if the addition 
would bring the result back into 
the representable range. 

Exponent overflow or exponent 
underflow during the multiplication 
is not recognized if the product in 
step 2 is set to a true zero. 

4. If instruction execution continues, 
the result element is stored at the 
location specified by general reg- 
ister 2. Then, if R t is zero, 8 is 
added to the contents of general 
registers 2 and 3; if Ri is zero and 
R 2 is not zero, 8 is also added to 
the contents of the R2 register. If 
Ri is not zero, the contents of the 
even and odd Rj registers are added 
to the contents of general registers 
2 and 3, respectively; if both Rx 
and R 2 are not zero, the contents of 
the even Rj register are also added 
to the contents of the R 2 register. 

5 . The contents of general register 1 
are decremented by one. If the 
result is zero, condition code is 
set, and instruction execution is 
completed. 

If the scalar factor in floating-point 
register is zero, it depends on the 
model whether an unnormalized element of 
operand B is recognized and condition 
code 2 is set in step 1, or whether 
instruction execution instead proceeds 
to step 2 and sets the product in step 2 
to zero. 

In step 2, the multiplication operation 
is the same as for the floating-point 
instruction MULTIPLY (MD) , and the addi- 
tion operation is the same as for ADD 
NORMALIZED (AD). Only the operand 
sources, the result target, and the han- 
dling of exception conditions differ. 



During the additions in step 4, carries 
out of bit position 8 for 24-bit 
I addressing and carries out of bit posi- 
| tion 1 for 31 -bit addressing are 
ignored; bit positions 0-7 or bit posi- 
tion 0, respectively, of the updated 
general registers are set to zeros . 

Floating-point register 0, the element - 
separation registers, and the elements 
in storage for operands B and C remain 
unchanged. 

Execution of the instruction is inter- 
ruptible by any interruption condition 
for which the CPU is enabled, other than 
an exponent-overflow, exponent- 
underflow, or significance exception. A 
unit of operation consists of one or 
more repetitions of the preceding five 
steps, with any interruption occurring 
at the end of step 5 or when the 
instruction is completed. 

When an interruption occurs during exe- 
cution and the interruption condition is 
not one that causes termination, general 
register 1 indicates the number of ele- 
ments remaining to be processed. The 
operand addresses have been updated to 
indicate the next set of elements to be 
processed, and the condition code is 
unpredictable . 

Access exceptions for operands may be 
recognized for storage locations other 
than the locations containing the 
current vector elements. For each 
operand, however, access exceptions are 
not recognized for more than one element 
beyond the current element or for 
element locations beyond the last 
element specified. 

I If any of the three tests made at the 
I start of instruction execution sets the 

I condition code, so that instruction exe- 

f cution is completed immediately, then no 
access exceptions are recognized for any 
operand, the change bits for operand A 
are unaffected, and no PER event for 
general -register alteration or storage 
alteration is indicated. 
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The storage location of operand A may 
coincide with the location of operand B 
or C, if the same first-element 
addresses and the same element sepa- 
rations are specified. If both condi- 
tions are not satisfied and partial 
overlap occurs between the location of 
operand A and the location of operand B 
or C, the contents of the location of 
operand A are undefined. Those contents 
are also undefined if R z = 2 or, for 
R 2 * 0, if: 

R 2 < 4, or 
R 2 = R a , or 
R 2 = Ri+1. 

Resulting Condition Code: 

All elements processed 

1 Exponent overflow, exponent under- 
flow, or significance loss 

2 Unnormalized scalar or operand-B 
element 

3 Element address or separation not 
multiple of 8 



Program Exceptions: 

• Access (fetch, operands B and C; 
store, operand A) 

• Operation (if the multiply-and-add 
facility is not installed) 

• Specification 



Programming Notes 

1. Unlike the scalar floating-point 
instructions MULTIPLY and ADD NOR- 
MALIZED, the MULTIPLY AND ADD 
instruction requires operands to be 
aligned on doubleword boundaries in 
storage. Moreover, unlike scalar 
MULTIPLY, it does not multiply 
unnormalized operands. The MULTIPLY 



AND ADD instruction also does not 
cause a program interruption when an 
arithmetic-exception condition is 
recognized for which the CPU is 
enabled. In all these cases, the 
instruction sets a nonzero condition 
code so that scalar floating-point 
instructions may be used to perform 
the arithmetic. 

If execution of the MULTIPLY AND ADD 
instruction sets a nonzero condition 
code, fewer than the specified 
number of vector elements were proc- 
essed. General register 1 contains 
the number of elements remaining to 
be processed in each vector. The 
general registers containing 
addresses designate the set of ele- 
ments which caused the exceptional 
condition. The program may then 
execute scalar floating-point 
instructions in an attempt to 
process these elements. 

If, during the execution of these 
scalar instructions, an exponent- 
overflow, exponent -underflow, or 
significance exception is recog- 
nized, and the corresponding program 
interruption is allowed, the inter- 
ruption may invoke standard fixup 
routines for these causes. If an 
unnormalized operand is encountered 
or a vector in storage is unaligned, 
the scalar instructions can process 
the elements . 

The following example in assembler 
language illustrates the type of 
programming that is recommended. By 
adding linkage and initialization 
instructions, a library subroutine 
is created that may be called from a 
high-level language, such as 
FORTRAN. 
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LI MADS 4,6 MULTIPLY AND ADD 

BC 8,L2 Test condition code 

* Do following instructions if CC not 

LD 2,0(0,3) Load element of B 

MDR 2,0 Multiply by S 

AD 2,0(0,6) Add element of C 

STD 2,0(0,2) Store element of A 

LA 2,0(4,2) Update address of A 

LA 3,0(5,3) Update address of B 

LA 6,0(4,6) Update address of C 

BCT 1,L1 Branch if not done 

* End of equivalent instructions 
L2 EQU * Continue 



In this example, the Rj and R 2 
fields of the MULTIPLY AND ADD 
instruction designate general regis- 
ters 4 and 6 . Thus , the number of 
elements is in general register 1; 
the addresses for operands A, B, and 
C are in general registers 2, 3, and 
6; and the element separations are 
in general registers 4 (for operands 
A and C) and 5 (for operand B). The 
MULTIPLY AND ADD instruction is fol- 
lowed by a BRANCH ON CONDITION 
instruction and by a loop containing 
equivalent scalar floating-point 
instructions, which are executed 
only when a nonzero condition code 
occurs. The loop uses floating- 
point register 2 as a working reg- 
ister. 

Note that this loop performs the 
same operations as one iteration of 
the preceding MULTIPLY AND ADD 
ins t ruct ion , except : 

• The storage operands may be una- 
ligned, and the element sepa- 
rations need not be a multiple 
of 8. 

• The multiplication operands may 
be in unnormalized form. 

• Exponent overflow causes an 
interruption; and exponent 
underflow and significance loss 
cause an interruption, if per- 



mitted by the program mask in 
the PSW. 

• A floating-point working reg- 
ister is needed. 

• The loop needs no initial test 
for a zero number of elements in 
general register 1, because the 
MULTIPLY AND ADD instruction 
sets condition code for this 
case. 

For the special instances in which 
vectors C and A coincide or in 
which successive vector elements 
are in successive doublewords, the 
routine may be simplified accord- 
ingly. 

3. Any two or all three vectors may 
coincide. The vectors coincide if 
both their addresses and their 
element separations are the same. 
Partial overlap in the storage areas 
for B and C may occur if neither 
vector overlaps with the result 
vector A. Partial overlap with the 
result vector, however, has unpre- 
dictable effects, which may differ 
from one model to another or from 
one execution to another. No check 
is made for partial overlap with the 
result vector. 

If vector C coincides with vector B, 
the same address should be loaded 
into general register 3 and into the 
R 2 register. The R 2 field should 
not designate general register 3; 
otherwise, it is unpredictable 
whether register 3 is updated once 
or twice for each set of elements 
processed. Similarly, the R 2 field 
should not designate general reg- 
ister 2 if vector C coincides with 
vector A. Thus, an R^ field con- 
taining 2 may or may not have the 
same effect as an R 2 field con- 
taining 0. 

4. Register and storage contents remain 
unchanged when instruction execution 
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ends before any elements have been 
processed successfully. 

See the section "interruptible 
Instructions" in Chapter 5, "Program 
Execution," of the appropriate Prin- 
ciples of Operation publication for 
more information concerning inter- 
ruptible instructions. Also, see 
the programming notes at the end of 
the section "Program-Event 
Recording" in Chapter 4, "Control," 
of that publication regarding redun- 
dant PER events that may occur when 
an interruptible instruction is 
resumed after an interruption. 

Special precautions must be taken if 
MULTIPLY AND ADD is made the target 
of EXECUTE. See the programming 
note concerning interruptible 
instructions under EXECUTE in 
Chapter 7, "General Instructions," 
of the appropriate Principles of 
Operation publication. 



The MULTIPLY AND ADD instruction is 
not equivalent to the loop con- 
taining scalar floating-point 
instructions in the above example if 
accessing an operand in storage 
causes an addressing or protection 
exception. If such an exception is 
recognized, execution of the MUL- 
TIPLY AND ADD instruction may be 
terminated such that the results 
left in the storage areas and in the 
general registers are unpredictable, 
and the condition code may not indi- 
cate an exceptional arithmetic con- 
dition for the elements being 
processed. Therefore, the program 
must not rely on the results and 
attempt to resume execution after an 
addressing or protection exception. 
(See the section "Termination" in 
Chapter 5, "Program Execution," of 
the appropriate Principles of Opera- 
tion publication.) 
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SQUARE-ROOT FACILITY 



The square-root facility consists of the 
SQUARE ROOT instruction and the square- 
root exception. The instruction, which 
extracts the square root of a floating- 
point operand in either the long or 
short format, is four bytes long and 
uses the RRE instruction format. The 
source operand resides in a floating- 
point register, and the result is placed 
in a floating-point register. When a 
positive, nonzero source operand is 
encountered which is unnormalized, it is 
first normalized at the start of the 
operation, but without changing the con- 
tents of the source-operand location. 
Nonzero results are always normalized. 
A zero result is made a true zero. 

The numeric result of the square-root 
operation on a valid operand is deter- 
mined by the rules of arithmetic. If 
the result can be represented exactly in 
the specified floating-point format, the 
exact result is produced. If the result 
cannot be represented exactly, it is 
rounded to the nearest number that is 
representable in the specified floating- 
point format. 

An operand that is less than zero is 
invalid. If the SQUARE ROOT instruction 
is executed with an invalid operand, a 
square-root exception is recognized. 

Square-Root Exception 

A square-root exception is recognized 
when the second operand of SQUARE ROOT 
is less than zero. 



INSTRUCTION DESCRIPTION 

The inclusion of instruction mnemonics 
in the description does not necessarily 
imply that the mnemonics are recognized 
by current assembler programs. A pro- 
grammer wishing to use this instruction 
in assembler- language programs with an 
assembler which does not recognize the 
mnemonics may use appropriate macros to 
assist with the translation of the 
instruction into machine language. 

SQUARE ROOT Instruction 



SQDR 



Ri,R 2 [RRE, Long Operands] 



I 'B244' | IIIIIIH I Ri I R 2 I 

I l I l l 

16 24 28 31 

SQER RijR2 [RRE, Short Operands] 



'B245' 1 1 1 II I III | R t | R 2 | 

I l I l 



16 



24 28 31 



The normalized and rounded square root 
of the second operand is placed in the 
first-operand location. 

When the fraction of the second operand 
is zero, the sign and characteristic of 
the second operand are ignored, and the 
operation is completed by placing a true 
zero in the first-operand location. 



The operation is suppressed. 

The instruction- length code is 2. 



When the second operand is less than 
zero, a square-root exception is recog- 
nized. 



The square-root exception is indicated 
by a program-interruption code of 00 ID 
hex (or 009D hex if a concurrent PER 
event is indicated) . 



When the second operand is normalized 
and greater than zero, the character- 
istic, fraction, and sign of the result 
are produced as follows: 
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• The result characteristic is one- 
half of the sum of the operand char- 
acteristic and either 64, if the 
operand characteristic is even, or 
65, if it is odd. 

• If the operand characteristic is 
odd, the operand fraction is shifted 
right one digit position, the right- 
most digit entering the guard-digit 
position. 

• An intermediate-result fraction is 
produced by computing without 
rounding the square root of the 
operand fraction, after any right 
shift as described. The 
intermediate-result fraction con- 
sists of the 15 most significant 
hexadecimal digits of the square- 
root result in the long format, or 
seven in the short format, where 
both formats include a guard digit 
on the right . 

• A one is added to the leftmost bit 
of the guard digit of the interme- 
diate result, any carry is propa- 
gated to the left, and the guard 
digit is dropped to produce the 
result fraction. 

• The result sign is made plus. 

When the second operand is unnormalized 
and greater than zero, the operand is 
first normalized. The operation then 
proceeds as for normalized operands. 

The R! and R 2 fields must designate reg- 
ister 0, 2, 4, or 6. Otherwise, a spec- 
ification exception is recognized. 



The code remains 



Condition Code : 
unchanged . 

Program Exceptions: 



Operation (if the square-root 
facility is not installed) 
Specification 
Square root 



Programming Notes 

1. The use of the SQUARE ROOT instruc- 
tion with short operands (SQER) is 
illustrated by the examples in 
Figure 4. 



1 

| Operand 


I 
Decimal | Result 


1 
Decimal | 


1 (hex) 

1 


Value | (hex) 

1 


Value | 


1 

|42 190000 


1 
25.0 |41 500000 


5.0 | 


|40 400000 


0.250 |40 800000 


0.50 | 


| 40 800000 


0.50 [40 B504F3 


0.7071. . . | 


|41 800000 

i 


8.0 | 41 2D413D 

i 


2.8284. . . | 

i 



Figure 4. Square-Root Examples 



2. The result fraction is correctly 
normalized without any further left 
or right shifts of the 
intermediate-result fraction and 
without any further exponent adjust- 
ment. Rounding cannot cause a carry 
out of the leftmost digit . 

3. Although a characteristic greater 
than 127 or less than zero may tem- 
porarily be generated during ( the 
operation, the result characteristic 
is always within the representable 
range, and no exponent overflow or 
underflow occurs. 

Specifically, the smallest nonzero 
operand in the long format consists 
of a one bit, preceded on the left 
by 63 zeros . This operand is an 
unnormalized number with a value of 
16" 78 , and its square root is 16" 39 . 
The normalized representation of 
this result has a characteristic of 
26 (decimal) . Similarly, the square 
root of the largest representable 

operand has a characteristic of % 

(decimal). The instruction, there- 
fore, cannot produce a nonzero 
result with a characteristic outside 
the range of 26 to 96. 
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MATHEMATICAL-FUNCTION FACILITIES 



The mathematical-function facilities 
include instructions which perform some 
of the mathematical operations that are 
commonly included as basic functions in 
higher-level programming languages. The 
instructions are ARCTANGENT, COMMON LOG- 
ARITHM, COSINE, EXPONENTIAL, NATURAL 
LOGARITHM, RAISE TO POWER, and SINE. 
They operate on floating-point numbers 
in either the long or short format . 

Each instruction is four bytes long and 
uses the RRE format . The source operand 
is in a floating-point register, and the 
result is placed in a floating-point 
register. When a nonzero source operand 
is encountered which is unnormalized, it 
is first normalized at the start of the 
operation, but without changing the con- 
tents of the source-operand location. 
Nonzero results are always normalized. 
A zero result is made a true zero. 

When a mathematical -function instruction 
performs a floating-point operation in 
the short format, the rightmost 32 bits 
of the source register are ignored as 
the operand is fetched. When the 
instruction places a short -format result 
in the target register, the rightmost 32 
bit positions of the register remain 
unchanged . 

The specific result produced by one of 
these instructions for a particular 
source operand or operand pair either 
equals or approximates the mathemat- 
ically exact result. When it is neces- 
sary to approximate the exact result, 
the instruction returns one of the two 
normalized floating-point numbers which 
are the nearest neighbors of the infi- 
nitely precise result; which of the two 
numbers is returned depends on the algo- 
rithm employed in the implementation. 
Thus, the maximum error is less than one 
unit in the last place of the normalized 
result. 



Arithmetic exceptions such as exponent 
overflow, operands outside the range 
implemented for a particular function, 
and other exceptional conditions do not 
cause a program interruption but are 
indicated by the condition code. If the 
operation is performed successfully, the 
instruction sets condition code 0; oth- 
erwise, the instruction a nonzero condi- 
tion code, as described for each 
instruction, and the result location 
remains unchanged. 

Different models may implement a dif- 
ferent selection of these instructions, 
or none at all. The operand range 
covered for a given instruction may not 
always be the same from one model to 
another. When an instruction operand is 
outside the operand range implemented by 
the model on which the instruction is 
executed, condition code 3 is set. 
Refer to the functional-characteristics 
manual of a particular model for the 
selection of instructions and their 
operand ranges . 

To avoid model -dependent operation, a 
program using one of the mathematical - 
function instructions should be accompa- 
nied by a subroutine which performs the 
entire function without relying on that 
instruction. Each use of the instruc- 
tion should always be followed by a con- 
ditional branch to test the condition 
code. If the code is not zero, indi- 
cating an exceptional condition such as 
operands that are outside the imple- 
mented range, a branch to the subroutine 
may then be used to complete the opera- 
tion or to handle the exception appro- 
priately. 



Programming Notes 

1. The instructions evaluate the fol- 
lowing mathematical functions: 
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r 
r 

r 
r 
r 



arctan(x) 
log(x) 

cos(x) 
exp(x) 
ln(x) 



r = y**x 
r = sin(x) 



for ARCTANGENT, 

for COMMON 

LOGARITHM, 

for COSINE, 

for EXPONENTIAL, 

for NATURAL 

LOGARITHM, 

for RAISE TO POWER, 

for SINE, 



where operand x is in the register 
designated as the second-operand 
location, operand y is in the reg- 
ister designated as the first- 
operand location, and ** indicates 
exponentiation. The result, r, is 
placed in the first-operand 
location, replacing the operand y, 
if any. 

2. NATURAL LOGARITHM and EXPONENTIAL 
perform inverse operations; that is, 
within operand ranges appropriate to 
each function, they approximate the 
relations : 

x = ln(exp(x)) and x = exp(ln(x)) 

Similarly, COMMON LOGARITHM and 
RAISE TO POWER with a base (first 
operand) of 10 perform inverse oper- 
ations which approximate the 
relations : 

x = log(10**x) and x = 10**(log(x)) 

3. The result of using RAISE TO POWER 
to evaluate y**0.5 for a positive y 
may differ slightly from the result 
of applying SQUARE ROOT to the same 
operand in the same floating-point 
format. Likewise, the result of 
using RAISE TO POWER to evaluate 
y**(±n), where y is positive and n 
is an integer, may not be the same 
as the result obtained by using MUL- 
TIPLY and DIVIDE instructions. 



Mathematical Constants 

In the instruction descriptions, e 
refers to the mathematical constant with 
the approximate value 2.718281828... and 
pi refers to the mathematical constant 
with the approximate value 
3.141592653... 



INSTRUCTION DESCRIPTIONS 
ARCTANGENT Instruction 

For long operands: 



| 'B294' | //////// 1 Ri | R 2 | 

l l l l I 



16 



For short operands: 



24 



28 31 



'B295' | //////// 1 Ri I R 2 I 

I I I I 



16 



24 



28 31 



The arctangent of the second operand is 
placed in the first -operand location. 
The operand and the result are 
floating-point numbers in the same 
format . 

The result is in radians, normalized, 
and smaller in magnitude than pi/2. A 
nonzero result has the same sign as the 
operand. If the operand has a zero 
fraction, the result is a true zero. 

If the normalized result would have a 
characteristic that is less than zero 
(exponent underflow), condition code 1 
is set; the setting of the exponent - 

underflow mask bit in the P5W has no 

effect. If the second operand is 
outside the operand range implemented by 
the model, condition code 3 is set. In 
both cases where a nonzero condition 
code is set, the first -operand location 
remains unchanged. 
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If the operation is completed normally, 
condition code is set. 

The Rj and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 

Valid operation 

1 Exponent underflow 
2 

3 Invalid operation 



Program Exceptions: 

• Operation (if the arctangent 
facility is not installed) 

• Specification 



Programming Note 

Exponent underflow can occur only when a 
nonzero second operand is unnormalized 
and so small that the result after nor- 
malization would have a characteristic 
less than zero. 

COMMON LOGARITHM Instruction 



The common logarithm (the logarithm to 
the base 10) of the second operand 
placed in the first-operand location. 
The operand and the result are 
floating-point numbers in the same 
format . 

The result is normalized. If the 
operand is plus one, the result is a 
true zero. 

If the second operand has a zero frac- 
tion or is negative, or if the second 
operand is outside the operand range 
implemented by the model, condition code 
3 is set, and the first -operand location 
remains unchanged. Otherwise, the oper- 
ation is completed normally, and condi- 
tion code is set. 

The Ri and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 




1 
2 

I 3 



Valid operation 
Invalid operation 



For long operands : 



'B284' | II HUH | Ri | R 2 I 

I I i l 



16 



24 28 31 



Program Exceptions: 

• Operation (if the common- logarithm 
facility is not installed) 

• Specification 



For short operands : 



'B285' | III///II | R x | R 2 

l I l 



16 



24 28 31 



Programming Note 

Exponent underflow cannot occur, because 
the result is either zero or suffi- 
ciently greater than zero in magnitude 
to be representable as a normalized 
floating-point number. 
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COSINE Instruction 

For long operands : 



'B28A' | //////// 1 Ri | R 2 | 

1 l l l 



24 



28 31 



16 



For short operands : 



'B28B' | //////// I Ri I R 2 

l I 1 



16 



24 



28 31 



The cosine of the second operand in 
radians is placed in the first -operand 
location. The operand and the result 
are floating-point numbers in the same 
format . 

The result is normalized and never 
greater in magnitude than one. If the 
operand has a zero fraction, the result 
is exactly plus one. 

If the second operand is outside the 
operand range implemented by the model, 
or if the absolute value of the second 
operand is not less than pi multiplied 
by 2 50 (long) or pi multiplied by 2 1 * 
(short), then the first-operand location 
remains unchanged, and condition code 3 
is set. Otherwise, the operation is 
completed normally, and condition code 
is set. 

The Rj and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 



Resulting Condition Code: 
Valid operation 



/ 

I 1 

I 2 

| 3 Invalid operation 



Program Exceptions: 

• Operation (if the sine-cosine 
facility is not installed) 

• Specification 



Programming Note 

Mathematically, the cosine of pi/ 2, or 
of any odd multiple of pi/2, is zero. 
In practice, the floating-point operands 
nearest to those values are sufficiently 
different from a multiple of pi/2 that 
the result cannot be zero or produce 
exponent underflow. 

EXPONENTIAL Instruction 

For long operands : 



*B282' I////////I Ri | R 2 | 

1 I I I 



16 



For short operands : 



24 



28 31 



'B283 1 



I II II I II 1 1 Ri 
j i 



16 



24 



28 31 



The result of raising the mathematical 
constant e to the power of the second 
operand is placed in the first -operand 
location. The operand and the result 
are floating-point numbers in the same 
format . 

The result is normalized. If the 
operand has a zero fraction, the result 
is exactly plus one. 

If the normalized result would have a 

characteristic that is less than zero 
(exponent underflow), condition code 1 
is set; the setting of the exponent- 
underflow mask bit in the PSW has no 
effect. If the normalized result would 
have a characteristic that is greater 
than 127 (exponent overflow) , condition 
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code 2 is set. If the second operand is 
outside the operand range implemented by 
the model, condition code 3 is set. In 
all three cases where a nonzero condi- 
tion code is set, the first-operand 
location remains unchanged. 

If the operation is completed normally, 
condition code is set . 

The R a and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 

Valid operation 

1 Exponent underflow 

2 Exponent overflow 

3 Invalid operation 



Program Exceptions: 

• Operation (if the exponential 
facility is not installed) 

• Specification 

NATURAL LOGARITHM Instruction 
For long operands : 



The natural logarithm (the logarithm to 
the base e) of the second operand is 
placed in the first-operand location. 
The operand and the result are 
floating-point numbers in the same 
format . 

The result is normalized. If the 
operand is plus one, the result is a 
true zero. 

If the second operand has a zero frac- 
tion or is negative, or if the second 
operand is outside the operand range 
implemented by the model, the first - 
operand location remains unchanged, and 
condition code 3 is set. Otherwise, the 
operation is completed normally, and 
condition code is set. 

The Ri and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 




1 
2 

I 3 



Valid operation 
Invalid operation 



'B280* | II I II Hi | Ri I R 2 | 

I I i I 



16 



24 



28 31 



Program Exceptions: 

• Operation (if the natural -logarithm 
facility is not installed) 

• Specification 



For short operands : 



'B281' I////////I Ri I R 2 I 
I i i I 



16 



24 



28 31 



Programming Note 

Exponent underflow cannot occur, because 
the result is either zero or suffi- 
ciently greater than zero in magnitude 
to be representable as a normalized 
floating-point number. 
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RAISE TO POWER Instruction 

For long operands : 



'B286' I//////// | Ri | R 2 | 

1 I i i 



16 



For short operands : 



24 



28 31 



'B287' 



j i 



1 1 

Ri I R 2 I 
i i 



16 



24 



28 31 



The first operand is raised to the power 
of the second operand, and the result is 
placed in the first-operand location. 
The operands and the result are 
floating-point numbers in the same 
format . 

The result is normalized. If the first 
operand has a zero fraction and the 
second operand is greater than zero, the 
result is a true zero. If the first 
operand has a nonzero fraction and the 
second operand has a zero fraction, the 
result is exactly plus one. 

For the operation to be valid, (1) the 
first operand must be greater than zero, 
or (2) if the first operand has a zero 
fraction, the second operand must be 
greater than zero, or (3) if the first 
operand is less than zero, the second 
operand must have a zero fraction. 

If the normalized result would have a 
characteristic that is less than zero 
(exponent underflow), condition code 1 
is set; the setting of the exponent - 



underflow mask bit in the PSW has no 
effect. If the normalized result would 
have a characteristic that is greater 
than 127 (exponent overflow) , condition 
code 2 is set. If the operation is 
invalid, or if either operand is outside 
the operand range implemented by the 
model, condition code 3 is set. In all 
three cases where a nonzero condition 
code is set, the first-operand location 
remains unchanged. 

If the operation is completed normally, 
condition code is set. 

The Rj and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 

Valid operation 

1 Exponent underflow 

2 Exponent overflow 

3 Invalid operation 



Program Exceptions: 

• Operation (if the raise-to-power 
facility is not installed) 

• Specification 



Programming Note 

Whether the operation is valid depends 
on the signs and on the zero or nonzero 
values of the two operands . Figure 5 
summarizes the conditions under which 
the operation is invalid, causing condi- 
tion code 3 to be set, or valid, causing 
condition code to be set except when 
exponent underflow (code 1) or exponent 
overflow (code 2) occurs. 
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1 1 1 
| | Second | 
| First | Operand j 


1 |< o|= o|> oj 


1 < 1 I 1 N | I | 


1 =0 | I | I | Z | 


| > | V | N | V | 

II! 


| Explanation : \ 

| I Invalid | 
| N Valid, result is one. | 
| V Valid | 
| Z Valid, result is zero.| 



Figure 5. Operand Validity 



I operand has a zero fraction, the result 

| is a true zero. 

| If the normalized result would have a 

| characteristic that is less than zero 

| (exponent underflow), condition code 1 

| is set; the setting of the exponent - 

| underflow mask bit in the PSW has no 

| effect. If the second operand is 

| outside the operand range implemented by 

| the model, or if the absolute value of 

| the second operand is not less than pi 

| multiplied by 2 50 (long) or pi multi- 

| plied by 2 1 ' (short), then condition 

| code 3 is set. In all those cases where 

| a nonzero condition code is set, the 

| first-operand location remains 

| unchanged . 

| If the operation is completed normally, 

| condition code is set. 



Not shown above is the possibility that 
condition code 3 is set because an oth- 
erwise valid operand is outside the 
range implemented by the model. 

SINE Instruction 

For long operands : 



The Ri and R 2 fields must designate reg- 
ister 0, 2, 4, or 6; otherwise, a spec- 
ification exception is recognized. 

Resulting Condition Code: 

Valid operation 

1 Exponent underflow 
2 

3 Invalid operation 



'B288* |//////// | Ri I R 2 | 

l l l I 



16 



For short operands : 



24 



28 31 



Program Exceptions: 

| • Operation (if the sine-cosine 
| facility is not installed) 
| • Specification 



•B289' I////////! Ri 



16 



24 28 31 



The sine of the second operand in 
radians is placed in the first-operand 
location. The operand and the result 
are floating-point numbers in the same 
format . 

The result is normalized and never 
greater in magnitude than one. If the 



Programming Note 

Mathematically, the sine of pi, or of 
any multiple of pi, is zero. In prac- 
tice, except for a zero operand, the 

floating-point operands nearest to those 

values are sufficiently different from a 
multiple of pi that the result cannot be 
zero or produce exponent underflow. 

Exponent underflow can occur only when a 
nonzero second operand is unnormalized 
and so small that the result after nor- 
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malization would have a characteristic 
less than zero. 



INSTRUCTION CHARACTERISTICS 

Figure 6 summarizes various character- 
istics of the mathematical-function 
instructions, including any restrictions 
on the operand and result range of each 
instruction, whether special conditions 
cause nonzero condition codes to be set, 



and the results produced by some 
instructions for certain unique oper- 
ands. All instructions set condition 
code when execution is completed 
normally, and all set condition code 3 
when an operand is outside the operand 
range implemented by the model. 



II III 
| | Range Restrictions on | Condition Code | | 


I 1 i 1 i i 1 1 
| Function | Operands | Result | 1 | 2 | 3 | Special Values | 


| ARCTANGENT j ||r| < pi/2\'E.\3 1 \ j jlf x=0, then r=0 | 
| COMMON LOGARITHM | x > j j | j x < | If x=l, then r=0 j 
JCOSINE j |x| < L ||r| < 1 | j j jlf x=0, then r=l j 
(EXPONENTIAL I | r > |EU |E0 j jlf x=0, then r=l j 
| NATURAL LOGARITHM! x>0 j j | x < |lf x=l, then r=0 


|RAISE TO POWER | (y>0)& any x| r > |EU |E0 | (y=0)&(x<0) | If (y=0)&(x>0) , r=0 | 
1 (y**x) |(y=0)&(x>0) j j | |(y<0)&(x*0)|lf(y*0)&(x=0), r=l j 
j j(y<0)&(x=0) j III 1 1 

JSINE j |x| < L ||r| < 1 j KU 1 1 j jlf x=0, then r=0 j 

i i i i i i 


\ Explanation: \ 

| & "and" | 
| EO Exponent overflow | 
| EU Exponent underflow | 
| EU 1 Exponent underflow can be caused only by a very small, unnormalized operand. | 
| L Limit: pi*2* ° for long operands, pi*2 l * for short operands | 
| r Result; placed in first-operand location | 
| r=0 Result is true zero. | 
| x Operand; obtained from second-operand location | 
| x=0 Second operand has zero fraction. | 
1 y Operand (only for some instructions); obtained from first-operand location | 
| y=0 First operand has zero fraction. I 



Figure 6. Instruction Characteristics 
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APPENDIX. INSTRUCTION SUMMARY 



1 1 ■ 1 1 

I 1 |0p | 
I Name | Characteristics |Code| 


1 1 i i l 1 1 
| ARCTANGENT (long) |RRE C ZT| SPJ j |B294| 
j ARCTANGENT (short) |RRE C ZT| SPJ j B295 j 
| COMMON LOGARITHM (long) |RRE C LCJ SPJ j B284| 
j COMMON LOGARITHM (short) |RRE C LC| SP j j | B285 | 
COSINE (long) |RRE C SN SP | | |B28A| 
COSINE (short) |RRE C SN | SP | | B28B j 


EXPONENTIAL (long) |RRE C EP| SP J | | B282 j 
j EXPONENTIAL (short) RRE C EP j SP j |B283| 
j MULTIPLY AND ADD |RRE C MA| A SP| II | R ST|B243| 
| NATURAL LOGARITHM (long) |RRE C LN| SP j | j B280 j 
JNATURAL LOGARITHM (short) |RRE C LN| SPJ |B28l| 


| RAISE TO POWER (long) |RRE C RP| SP | | B286| 
j RAISE TO POWER (short) |RRE C RP| SPJ j | B287 j 
|SINE (long) JRRE C SN SP | j j B288 | 
JSINE (short) |RRE C SN SP j | B289 | 
j SQUARE ROOT (long) |RRE QRJ SP SQ j |B244| 
| SQUARE ROOT (short) |RRE QR SP SQ j |B245| 

11111 


| Explanation : I 

| A Access exceptions for logical addresses | 
| C Condition code is set. I 
| EP Exponential facility I 
| II Interruptible instruction I 
| LC Common- logarithm facility I 
| LN Natural -logarithm facility I 
| MA Multiply -and -add facility I 
| QR Square -root facility I 
j R PER general-register-alteration event | 
| RP Raise-to-power facility I 
| RRE RRE instruction format I 
| SN Sine -cosine facility 1 
| SP Specification exception 1 
| SQ Square-root exception 1 
1 ST PER storage-alteration event I 
| ZT Arctangent facility 1 



Figure 7. Summary of Mathematical-Assist Instructions 
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access exceptions for MULTIPLY AND ADD 

6 
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boundary alignment of vector elements 1 
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change bit for MULTIPLY AND ADD 6 

COMMON LOGARITHM instructions 15 

condition code 

for mathematical functions 13 
for MULTIPLY AND ADD 2 

contiguous vector elements 1 

COSINE instructions 16 
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interruptible instruction 2 



logarithm 
common 15 
natural 17 



e (mathematical constant) 14 

error bounds 
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for SQUARE ROOT 11 

exceptions 
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square -root 11 
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MADS (MULTIPLY AND ADD) instruction 

mathematical constants e and pi 14 

matrix computations 1 

mode of operation 3 

MULTIPLY AND ADD (MADS) instruction 
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NATURAL LOGARITHM instructions 17 
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overflow 
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for RAISE TO POWER 18 



PER (program-event recording) for MUL- 
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pi (mathematical constant) 14 



underflow 

for ARCTANGENT 14 

for EXPONENTIAL 16 

for MULTIPLY AND ADD 2 

for RAISE TO POWER 18 

for SINE 19 
unit of operation for MULTIPLY AND ADD 

6 
unnormalized operands 

for mathematical -function 
instructions 13 

for MULTIPLY AND ADD 1 

for SQUARE ROOT 11 
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vector iii 



scalar iii 

SINE instructions 19 

SQDR (SQUARE ROOT) instruction 11 

SQER (SQUARE ROOT) instruction 11 

SQUARE ROOT (SQDR, SQER) instructions 

square-root exception 11 
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